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Investigation  of  the  Laser  Material  Interaction  Regimes 
During  Thermal  and  Compositional  Superposition 
In  Laser  Surface  Improvement 


BACKGROUND 


During  Laser  Surface  Alloying  (LS  A),  additives  are  applied  to  a  base  material  in  such  a  way 
that  upon  melting  a  thin  region  of  the  surface,  significant  alloying  between  the  additive  and 
the  base  material  occurs.  This  is  demonstrated 
in  Figure  1  which  is  a  schematic  of  a  LISF*^  (La¬ 
ser  Induced  Surface  Improvement)  surface.  In 
this  case,  the  additive  is  applied  with  a  surface 
sluny  prior  to  melting.  The  laser  beam  is  then 
passed  across  the  surface,  allo3dng  the  additive 
into  the  substrate.  To  cover  large  surfaces,  the 
laser  beam  is  indexed  in  sequential  overlaps.  The 
goal  of  LS  A  is  to  obtain  properties  that  both  dif¬ 
fer  fi'om  the  substrate,  and  provide  some  form  of 
desired  protection  (wear,  corrosion,  etc.).  This 
enables  the  user/designer  to  select  a  base  mate¬ 
rial  for  its  desirable  properties  without  having  to 
compromise  those  properties  in  order  to  survive 
a  harsh  environment. 

Since  LSA  involves  the  spatial  and  temporal  input  of  heat  into  the  base  material  in  order  to 
affect  the  alloying,  both  the  alloying  process  itself  and  the  influence  of  the  heat  upon  the 
surface  and  base  material  must  be  considered.  A  schematic  for  a  typical  material  is  seen  in 

Figure  2  and  shows  that  in  addition  to  the  sur¬ 
face  alloying  and  depending  upon  the  type  of 
base  material  (steel,  aluminum  alloy,  etc.),  nu¬ 
merous  metallurgical  events  can  occur  beneath 
the  alloyed  region.  Add  to  this  the  additional 
complexity  of  indexing  or  overlapping  succes¬ 
sive  beams  to  create  large  surfaces,  and  the 
challenge  of  understanding  the  LSA  process 
grows  to  tremendous  proportions.  This  re¬ 
search  seeks  to  address  some  aspects  of  that 
challenge. 

Figure  2.  A  schematic  of  the  regions  that 
occur  in  a  material  during  laser  surface. 
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Figure  1.  Schematic  ofLISI  Process. 
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OBJECTIVES 


The  objectives  of  the  research  are  to  1)  understand  the  reactions  occurring  within  the  in¬ 
dexed  and  surrounding  regions  as  a  result  of  the  thermal  and  compositional  cycling  experi¬ 
enced  during  the  formation  of  a  laser  surface  modified  layer,  and  2)  to  evaluate  the  process 
kinetics,  solidification  and  growth  morphologies  based  on  thermodynamic  and  heat-conduc¬ 
tion  considerations.  In  order  to  accomplish  these  objectives,  the  laser  surfacing  process  was 
divided  into  three  significant  categories.  These  are  1)  the  alloyed  surface  region  and  the 
phase  reactions  within  that  surface,  2)  the  grain  growth/solid  state  transformation  and  tem¬ 
pered  regions  and  the  transformations  within  those  regions,  and  3)  the  thermal  temporal 
history  occurring  during  processing  as  determined  both  ejq)erimentally  and  computationally. 

Processing  parameters  for  all  of  the  trials  are  given  in  Table  I.  Compositions  of  the  base 
alloys  used  are  given  in  Table  II.  In  all  cases,  a  Nd:YAG  fiber  optically  delivered  laser  was 
used. 


TABLE  I 

Processing  Conditions 


CHROMIUM  ADDITIONS 


1500 

4500 

0.25 

1.0 

0.50 

0.75 

1.0 

CHROMIUM,  NICKEL,  WC  AND  SiC  ADDITIONS 


Index  tmml 

2000 

1000 

2.5 

0.75  X  4 

1500 

2250 

3000 

CHROMIUM  AND  CrB2  ADDIHONS 


Sneed  tmm/sec) 

Index  (mm) 

1500 

2000 

3.5 

0.75x4 

1500 

TABLE  I  Cont. 


1000 

1.5 .2.5 

1250 

1. 5.2.5 

1500 

1.5 .2.5 

1750 

1.5 .2.5 

2000 

1.5. 2.5 

2250 

1.5 .2.5 

2500 

1.5 .2.5 

3000 

1.5 .2.5 

3500 

1.5 .2.5 

4000 

1.5 .2.5 

4500 

1. 5.2.5 

HEAT  TREATED  4340  PROCESSING 


Sneed  fmm/sec 


500 


1000 


Power  fwatts’ 


500 


750 


1000 


1500 
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THERMAL  PROCESSING  TRIALS 


Beam  Translation 
Velocity  (mm/min) 

Beam  Index  (mm) 

1000 

1.5,  2.5 

1250 

1.5,  2.5 

1500 

1.5, 2.5 

1750 

1.5,  2.5 

2000 

1.5,  2.5 

2250 

1.5,  2.5 

2500 

1.5,  2.5 

3000 

1.5,  2.5 

3500 

1.5,  2.5 

4000 

1.5,  2.5 

4500 

1.5,  2.5 

Beam  Size  (mm) 


1  x3 


TABLE  II 


Set 

Speed 

(mm/min) 

Layer 

Depth 

(pm) 

Average 

Hardness 

(Knoops) 

Average  Composition 
(rem  Fe) 

Total 

Carbides 

(wt%) 

Carbide/Iron/Silcon 
Phase  (1  ♦►4) 

Cr 

wt% 

Ni 

wt% 

Si(C) 

wt% 

W(C) 

wt% 

A 

1000 

462 

576  ±64 

10.2 

2.0 

7.5 

4.2 

11.7 

(Least)  1 

B 

1500 

315 

634  ±59 

12.4 

2.0 

8.7 

5.2 

13.9 

2 

C 

2250 

231 

644  ±97 

13.3 

2.0 

9.1 

6.3 

15.4 

3 

D 

3000 

246 

741  ±161 

17.9 

2.0 

12.2 

8.4 

20.6 

(Most)  4 

E  (4340P) 

1500 

590 

0.8 

1.8 

0 

0 

0 

F  (4340) 

270 

IL-6096 

L  The  Alloyed  Surface  Region 

A  systematic  approach  was  adopted  to  investigating  the  alloyed  surface  region  in  steels. 
This  began  with  the  addition  of  one  element,  chromium,  and  continued  with  additions  of 
nickel,  WC  and  SiC,  and  then  CrB^. 

Chromium  Additions 

A  series  of  sequential  laser  surfece  modification  experiments  was  completed  with  the  associ¬ 
ated  compositional,  microstructural  and  mechanical  property  determinations.  This  involved 
a  matrix  of  samples  which  were  processed  on  1 0 1 0  steel  with  chromium  additive  since  there 
is  a  large  amount  of  information  on  iron/chromium  alloys.  The  resultant  alloyed  regions 
ranged  in  depth  fi-om  100  microns  to  500  microns  and  with  total  widths  up  to  60  mm. 
Selected  track  sequences  on  each  sample  were  submitted  to  examination  at  ten  micron  inter¬ 
vals  to  produce  maps  of  composition,  phase  and  microhardness.  The  compositions  matched 
those  calculated  using  a  precursor  density  of  80%.  The  phases  were  predictable  by  using  the 
Schaeffler  diagram  [1].  This  indicates  that  for  LSA  using  single,  non-complex  elements, 
predictions  of  composition  and  phase  are  straight  forward. 

Surface  roughness  variation  was  dependant  on  processing  energy  while  the  layer  depth  and 
volume  for  each  alloy  was  found  to  be  dependent  not  only  on  input  power  and  energy  but  on 
the  indexing.  As  the  index  increased  (Le.,  the  distance  between  successive  passes  increased), 
the  depth  and  volume  increased.  This  is  due  to  the  fact  that  when  the  leiser  is  indexed,  it 
strikes  both  already-processed  surface  and  new  surface  slurry.  Since  an  increase  in  index 
assures  that  there  will  be  a  larger  area  of  fi-esh  slurry,  the  absorptivity  to  the  Nd:  YAG  laser 
beam  is  being  assisted  by  the  presence  of  the  slurry.  Other  possible  contributors,  such  as 
changes  in  conductivity  or  fluidity  due  to  the  chromixmi  additions,  can  be  eliminated  sinee 
there  is  no  indexing  effect  on  the  amount  of  chromium  dissolved  per  voliune. 
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Layer  hardness  was  also  found  to  be  related  to  indexing  as  well  as  input  power.  This  is 
illustrated  in  Figure  3  which  shows  the  region  of  maximum  hardness  in  the  heat  affected 
zone  in  the  low  energy  (23.33  joules/mm^)  sample  while  the  higher  energy  (93.3  joules/mm^) 
sample  has  maximum  hardness  within  the  layer  itself  This  effect  is  felt  to  reflect  the  strong 
influence  of  repetitive  temperature  passes  on  both  the  alloyed  region  and  the  HAZ  and  em¬ 
phasizes  the  necessity  for  further  study  of  the  superposition  of  overlapping  tracks. 
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Figure  3. 
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Figure  3.  Cont. 


Chromium.  Nickel.  WC  and  SiC  Additions 

Further  additions  were  made  to  the  surface  slurry  in  order  to  alter  the  composition  and 
therefore  phase  transformations.  The  addition  of  the  elements  chromium  and  nickel  enhance 
the  corrosion  resistance,  particularly  when  chromium  exceeds  12  wt%,  making  them  desir¬ 
able  additions  to  the  surface  of  a  steel  material.  They  also  dictate  predictable  phase  changes, 
with  chromium  being  a  ferrite  former  and  nickel  an  austenite  former,  therefore,  reducing  the 
hardening  due  to  martensitic  transformation.  The  phases  present  should  be  predictable  from 
the  Schaeffler  diagram  using  compositional  analysis  of  the  layers  and  the  overall  hardness 
can  then  be  predicted  using  weighted  hardness  values  for  the  individual  phases.  These  phase 
alterations  can  produce  a  significant  secondary  effect  in  the  surface  alloy  by  reversing  the 
residual  stresses  from  the  compressive  stresses  expected  in  a  martensitic  surface  to  the  ten¬ 
sile  stresses  in  an  austenite  surface  [2, 3].  Corrosion  protection  can  be  obtained  while  pro¬ 
ducing  a  variety  of  steel  phases  by  varying  the  composition  with  the  processing  parameters. 

Carbides  on  the  other  hand  increase  the  hardness  and  wear  resistance,  but  their  chemistry 
must  be  considered  during  processing  since  their  decomposition  can  alter  the  steel  phases  by 
supplying  either  carbon,  which  is  an  austenite  former,  the  individual  elements  of  tungsten 
and  silicon,  which  are  ferrite  formers,  or  new  carbides.  The  dissociation  of  the  carbide 
additives  is  a  strong  possibility  during  laser  processing  with  its  intense  energy  density.  For 
example,  although  WC  is  chemically  stable  rather  stable,  it  decomposes  into  its  components 
above  2600"C  [4],  a  temperature  well  within  the  regime  of  laser  surface  processing.  There¬ 
fore  the  WC  and  SiC  were  added  to  determine  if  the  surface  alloy  properties  could  be  pre¬ 
dicted  using  the  Schaeffler  diagram  and  hardness  values  for  the  phases  and  carbides. 

Steel  coupons  of  4340  steel  (6  x  76  x  154  mm)  were  used  as  substrates  for  these  experi¬ 
ments.  The  4340  was  selected  because  it  is  generally  the  standard  to  which  other  high 
strength  steels  are  compared  and  combines  hardenability  with  ductility,  toughness  and  strength. 
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In  addition,  it  has  good  strength  at  elevated  temperatures  and  is  readily  heat-treatable.  The 
large  surfaces  of  the  substrates  were  sandblasted  and  then  coated  with  a  slurry  comprised  of 
chromium  (44  wt%),  nickel  (22  wt%),  silicon  carbide  (3  wt%)  and  tungsten  carbide  (8  wt%) 
powder,  varnish  and  mineral  spirits.  Once  the  slurry  had  dried  at  room  temperature,  pro¬ 
cessing  proceeded  For  baseline  purposes,  an  uncoated  sample  of 4340  was  processed. 

Subsequent  to  processing,  the  sample  surfaces  were  examined  optically  for  signs  of  crack¬ 
ing.  Each  was  then  cut  into  pieces  for  more  extensive  analysis.  Energy  dispersive  spectros¬ 
copy  (EDS)  (standardless,  background  removed)  analysis  was  used  for  determination  of  the 
chromium,  nickel,  silicon  and  tungsten  content.  Using  a  second  set  of  samples,  fiirther 
anafysis  was  conducted  by  x-ray  diffraction  (CuKa,  1 .54  A°  wavelength  radiation),  first  on 
the  as-processed  surface,  then  on  the  same  with  100  microns  of  the  surface  removed.  Wear 
properties  of  were  ascertained  by  two  sets  of  wear  testing  (block-on-disc  apparatus  at  1000 
rpm  (280  mm/min)  and  4  lbs.  Normal  weight)  to  determine  weight  loss/minute  and  fi-iction 
coefficient  after  2, 10  and  20  minutes. 

A  compilation  of  the  characterization  results  is  given  in  Tables  11  and  III  (average  values). 
The  samples  (A,  B,  C,  and  D)  are  presented  in  order  of  increased  processing  speed. 
Comparison  data  for  4340  processed  (E)  and  as  received  (F)  is  also  given.  AH  of  the  sample 
surfaces  exhibited  some  evidence  of  surface  micro-cracking  (Figure  4)  indicative  of  residual 
stresses.  From  Table  II  it  can  be  seen  that  the  chromium  and  carbide  contents  vary  inversely 
to  the  depth  of  the  layers.  The  nickel  content  however,  remains  at  2%,  which  is  around  the 
limit  of  accuracy  of  the  analysis  technique.  The  cracking  increased  in  severity  fi’om  Set  A  to 
D  (Table  III),  with  D  containing  the  highest  amount  of  inclusion  induced  fl'acture  sites.  Set 
A  had  relatively  little  cracking  which  occurred  primarily  at  the  midpoint  of  the  beam  track. 
Sets  C  and  D  had  extensive  cracking  exhibiting  extremely  brittle  failure  and  penetrating  the 
depth  of  the  processed  region.  Using  calculations  of  the  chromium  and  nickel  equivalents 
and  a  calculated  value  of  the  amoimt  of  carbon  present  (assuming  dilution  fi-om  alloying),  the 
percentage  of  the  martensite,  austenite  and  ferrite  phases  in  each  sample  were  estimated 
fi-om  the  Schaeffler  diagram  and  are  given  in  Table  III.  This  was  contrary  to  what  was 


TABLE  m 


Set 

Cracking 
Severity 
(1-  4) 

Crack 

Location 

Predicted  Schaeffler 
Phases 

Observed  Phases 

Martensite 

vol% 

Ferrite 

vol% 

Austenite 

vol% 

A 

Overlap 

91 

0 

9 

Austenite/ferrite 

B 

Mid  track 

68 

0 

32 

Austenite/ferrite 

C 

3 

Both 

60 

0 

40 

Austenite/ferrite 

D 

(Most)  4 

Both/chunks 

4 

86 

Austenite/ferrite 

0 

Martensite/ferrite 
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Figure  4.  Surface  cracking  after  processing  at  2000  W  and  (a)  1000  mm/min,  sample  A  and  (b) 
2250  mm/min,  sample  C. 


observed  from  the  x-ray  diffraction  analysis  which  showed  that  A  contained  significant 
austenite  indicating  a  shift  toward  higher  carbon,  while  the  remaining  samples  had  noticeable 
ferrite  content.  Thus  the  predictions  using  the  Schaeffler  diagram  should  be  shifted  toward 
higher  chromium  equivalents  and  therefore  more  ferrite.  This  is  probably  due  to  the 
dissociation  of  the  carbides  since  indications  from  the  x-ray  of  any  remaining  undissociated 
SiC  and  WC  were  very  weak,  although  sample  D  evidenced  a  large  amount  of  SiC-5H.  In 
addition,  when  the  surfaces  were  removed,  the  remaining  layers  all  strongly  indicated  the 
presence  of  a  carbon-iron-silicon  (12,  79,  9  wt%)  phase  suggesting  a  dissociation  of  the 
incorporated  silicon  carbide  and  a  non-stoichiometric  composition  of  Si-C  compound.  The 
degree  of  occurrence  of  this  newly  formed  carbide  is  given  in  Table  11  and  indicates  the 
possibility  that  some  elemental  silicon,  which  is  a  ferrite  former,  was  present  in  the  alloyed 
layer.  Interestingly,  the  occurrence  of  the  C-Fe-Si  phase  increases  as  the  interaction  time 
(inverse  to  speed  of  processing)  decreased,  suggesting  a  non-equilibrium  phase. 

As  can  be  seen  in  Table  III,  the  average  microhardness  varies  inversely  with  layer  depth.  The 
standard  deviation  is  also  given,  so  an  appreciation  can  be  gained  for  the  extent  of  inhomo¬ 
geneity  within  the  sets,  which  is  also  reflected  in  the  microhardness  maps  shown  in  Figure  5 
(the  light  regions  indicate  high  hardness).  A  plot  of  average  microhardness  vs.  the  total 
silicon  plus  tungsten  percent  (Figure  6)  reveals  a  linear  relationship  which  when  curve-fit 
yields  the  following; 


Microhardness  =  373  +  17.9  (wt%  carbides) 

Using  this  relationship  between  hardness  and  weight  percent  silicon  plus  tungsten,  the  hardness 
of  the  silicon  plus  tungsten  additives  can  be  calculated  to  be  2100  Knoops,  an  acceptable 
value  for  carbides  based  on  the  literature  [5].  The  matrix  material  however  (zero  wt% 
Si+W)  has  a  Knoops  hardness  significantly  lower  than  the  hardness  of  sample  E,  the  processed 
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Wei^t  Percent  Total  Silicon  Plus  Tungsten 


Figure  6.  Plot  of  average  Knoops  microhardness  vs.  silicon  plus  tungsten  weight  percent. 
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matrix  material  with  no  aUoying  addition.  Sample  E’s  value  of  590  Knoops  is  in  keeping 
with  predictions  for  transformed  4340  with  approximately  100%  martensite.  This  suggests 
that  the  microhardness  of  the  resolidified  matrix  is  also  being  aflFected  by  the  alloying  additions 
and  the  formation  of  the  new  C-Fe-Si  phase  and  that  the  hardness  of  the  layer  is  dominated 
by  the  addition  of  the  carbides. 

Wear  test  results  are  given  in  Figure  7,  which  gives  grams  lost  per  minute  over  a  ten  minute 
period  for  the  base  4340(F)  as  well  as  the  processed  4340(E)  and  the  sample  surfaces.  In 
two  of  the  samples  (A  and  D),  the  addition  of  chromium,  nickel  and  the  carbides  increased 
the  wear  over  the  processed  and  unprocessed  base  material.  In  samples  B  and  C,  the  wear 
was  slightly  or  similar  to  the  base.  High  wear  would  be  expected  in  the  softer  sample  A,  and 
high  wear  would  also  be  ejqtected  in  the  harder  sample  D  which  was  seen  to  exhibit  signifi¬ 
cant  cracking  and  loss  of  chunks  of  material  during  testing.  This  is  not  expected  either  from 
the  Schaeffler  predictions  or  the  observed  phases  and  therefore  substantiates  the  conclusion 
that  there  are  significant  (and  sometimes  unknown)  phase  variations  during  laser  surface 
processing  when  a  variety  of  additives  are  employed. 


Figure  7.  Wear  slopes  of  processed  and  unprocessed  4340  samples  (block-on-disc  appa¬ 
ratus  at  280  m/min  and  4  lbs  normal  weight). 
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Chromium  and  CrB.  Additions 

Heat  treated  4340  (450  Knoops)  was  also  used  for  these  experiments.  Ten  micron  size  Cr 
and  CrB2  powders  were  mixed  in  ratios  from  lOCr/OOCrB^  to  90Cr/10CrB2  in  increments  of 
ten.  These  mixtures  were  then  dispersed  in  a  slurry,  sprayed  onto  the  steel  surface  to  a 
thickness  of  100  microns  and  processed.  The  samples  were  sectioned  for  compositional 
analysis,  microhardness,  surfece  roughness,  wear  testing  (block-on-disc  at  280  mm/min), 
abrasion  testing  (50  -  70  grit  A.F.S.  testing  sand  at  2000  revolution  increments,  wheel  diam¬ 
eter  8  inches)  and  x-ray  diffractometry  (CuK^  radiation).  Corrosion  properties  were  ascer¬ 
tained  by  corrosion  testing  in  a  Q-Fog  chamber  (ptt  6.5-7.2,  5  Avt%  NaCl,  35°C,  100% 
relative  humidity)  for  24  hours. 

Data  from  some  of  the  analyses  are  shown  in  Table  IV.  Figures  8  through  1 0  show  the  results 
of  the  microhardness,  wear  and  abrasion  tests  as  a  function  of  the  percent  of  CrB^  in  the 
slurry.  The  microhardness  increases  linearly  with  increasing  CrB^  irntU  reaching  the  highest 
percentage  sample.  A  linear  regression  curve-fit  of  the  data  gives  an  average  slope  of  16 
Knoops  per  wt%  CrB^  and  an  intercept  of  167  Knoops.  The  low  intercept  value  indicates 
that  the  matrix  is  significantly  softer  than  quenched  4340  steel  and  probably  represents  the 
steel  alloyed  with  chromium  to  produce  softer  ferrite.  The  effect  of  the  hardness  increase 
can  be  seen  in  the  wear  test  results,  which  show  a  related  decrease  in  wear  as  the  CrB^ 
addition  increased.  While  the  abrasion  tests  were  not  as  consistent,  they  still  show  the  same 
trend  of  a  decrease  in  wear  with  increasing  CrB^.  All  of  the  samples  placed  in  the  Q-Fog 
showed  some  corrosion,  particularly  on  the  edges  since  they  were  not  protected.  The  samples 
with  90  to  50  wt%  chromium  in  the  slurry  had  very  minimal  corrosion  indicating  that  there 
was  enough  chromium  present  in  the  iron  alloy  to  fi)rm  a  protective  coating.  Samples  with 
less  chromium  (and  therefijre  having  60  to  90  wt%  CrB2)  evidenced  extensive  corrosion. 
Results  from  the  x-ray  diffraction  patterns  indicated  the  formation  of  additional  borides,  the 
possibility  some  chromium  carbides  (Cr^Cj  and  Cr23Cg)  as  reported  in  previous  work  [6]  and 
the  absence  of  the  initial  boride  CrB2.  Borides  found  were:  CrB^,  CrB,  (Cr,  Fe)2  B  and 
possibly  FOjB.  Since  a  laser  surfrice  processed  area  often  solidifies  with  preferred  orienta¬ 
tion,  it  was  not  possible  to  quantitatively  compare  the  amoimts  of  the  phases.  However, 
CrB^  and  CrB  were  more  dominant  in  the  samples  with  higher  amounts  of  chromium  while 
Fe2B  was  more  dominant  in  the  samples  with  the  higher  amounts  of  added  CrB2.  This 
indicates  when  chromium  was  more  available  the  reactions  were  driven  towards  those  com¬ 
pounds  that  are  higher  in  chromium.  In  contrast,  when  there  was  minimal  or  no  chromium 
alloyed  with  the  CrB2,  then  the  reactions  were  driven  towards  those  with  less  chromium  and 
containing  iron.  An  additional  sample  was  processed  that  contained  only  CrB2  in  the  slurry. 
SRD  analysis  on  the  surface  revealed  primarily  CrB2,  substantiating  the  conclusion  that  the 
presence  of  additional  chromium  drives  the  reactions  toward  other  chromium  borides.  The 
iron  boride  phases  were  not  observed  in  the  absence  of  chromium.  A  san^le  processed  with 
only  chromium  in  the  slurry  produced  the  Fe-Cr  phase  and  some  evidence  of  chromium 
carbides  (Cr^Cj  and  Cr23Cg).  X-ray  diffraction  patterns  taken  after  the  wear  tests  showed  no 
significant  change  in  the  phases  present  with  the  possible  occurrence  of  some  Fe^B. 


14 


Speed  Hardness  Dry  Wear  Abrasion  Layer  Surface 

(mm/min)  (Knoops)  (gms/min)  (gms/rev)  Corrosion  ^^^jCross-Sectioi^^^^  wt  %Cr  wt  %Cr 

90Cr/10CrB2  1000  351  0.0009  2.81167E-05  LowRust  ok  22.6  67.0 


TABLE  IV 

Analysis  and  Test  Results 
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Figure  8.  Microhardness  as  a  function  of  the  added  CrB^for  two  processing  speeds. 


Figure  9.  Dry  wear  rates  as  a  function  of  the  added  CrB^for  two  laser  processing  speeds. 


Figure  10.  Abrasion  wear  rates  as  a  function  of  the  added  CrB^for  two  laser  processing  speeds. 


Figures  11  through  13  show  representative  microstructures  from  three  of  the  samples.  Fig¬ 
ure  11  contains  the  most  precursor  chromium  (90  wt%)  and  least  CrB^  (10  wt%).  The 
solidification  structure  is  uniformly  dendritic  with  the  microhardness  varying  by  no  more 
than  ten  Knoops  from  the  surface  to  the  interface.  Figure  12  shows  the  structure  of  a  sample 
processed  with  50  wt%  Cr  and  50  wt%  CrB^  in  the  precursor.  The  larger  white  particle  has 
a  hardness  of  1718  Knoops  indicating  that  it  might  be  a  partially  dissociated  CrB^  particle 
which  should  have  high  hardness.  The  rest  of  the  matrix  varies  from  1000  to  2000  Knoops. 
The  dark  band  at  the  interface  has  a  hardness  of 375  Knoops  and  is  due  to  the  initial  solidi¬ 
fication  occurring  epitaxially  on  the  substrate.  In  regions  where  the  solid  particles  occur,  the 
white  trunks  are  seen  to  nucleate  from  those  particles.  This  is  shown  in  a  higher  magnified 
photograph  in  Figure  14.  Figure  13  shows  a  view  of  the  sample  with  20  wt%  chromium  and 
80  wt%  CrBj  in  the  precursor.  The  boride  columns  are  considerably  more  dense  than  in 
Figure  13  as  is  evidenced  by  the  microhardness  measurements  of  1640  Knoops  in  those 
regions. 


Figure  11.  Microstructure  of  laser 
processed  layer  (1500  W 1000  mm/min) 
with  90%  Cr/10%  CrB ^precursor 
addition. 


Figures  12.  Microstructure  of  laser 
processed  layer  (1500  W 1500  mm/ 
min)  with  50%  Cr/50%  CrB^  precursor 
addition. 
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Figure  13.  Microstructure  of  laser 
processed  Ic^er  (1500  W 1500  mm/ 
min)  with  20%  Cr/80%  CrB^precursor 
addition. 


(a) 

Figure  14a  and  b.  Microstructure  of 
laser  processed  layer  (1500  W 1500 
mm/min)  with  50%  Cr/50%  CrB^ 
precursor  additions  showing  growth 
of  columnar  structure  from  large 


n.  Solid  State  Transformation  Regions 


4340  Processin2  Trials 

A  significant  focus  of  the  effort  is  the  delineation  of  mechanical  properties  based  on  process¬ 
ing  conditions.  Therefore,  a  sequence  of  experiments  without  slurry  additions  was  con¬ 
ducted  in  which  coupons  of  6  mm  thick  4340  plate  were  processed  using  a  Nd:  YAG  laser 
beam  with  dimensions  of  3  mm  x  1  mm  and  a  power  of 500  Watts.  As  discussed  earlier,  4340 
steel  is  a  high  strength  structural  material  that  is  heat  treated  to  achieve  desired  properties. 
As  a  result,  when  this  material  is  laser  surface  alloyed,  there  are  considerable  metallurgical 
phenomena  that  can  occur  beneath  the  molten  region  [7,  8].  The  beam  translation  rate  and 
index  were  the  primary  processing  variables  (see  Table  I).  Two  sets  of  trials,  each  with  a 

different  beam  translation  index  be¬ 
tween  consecutive  tracks,  were  com¬ 
pleted.  Variations  in  index  yielded  dif¬ 
ferent  amounts  of  overlap  between  the 
adjacent  processed  regions  but  within 
the  beam  width.  The  resultant  sur¬ 
face  modified  regions  ranged  in  maxi¬ 
mum  depth  from  240  microns  to  300 
microns,  with  total  widths  up  to  35  mm 
(or  ten  laser  tracks).  Selected  track 
sequences  on  each  sample  were  sub¬ 
mitted  to  examination  at  100  micron 
intervals,  as  shown  in  Figure  15,  to 
produce  maps  of  microhardness. 

Figures  16  and  17  are  microhardness 
contour  maps  for  beam  translation 
rates  of  1000  mm/min  and  4500  mm/ 
min,  respectively,  for  the  beam  index 
of  1 .5  mm.  Although  the  laser  energy 
does  not  initiate  melting,  the 
microhardness  results  for  the  lower 
translation  rate  cases  indicate  that  a  transition  from  ferrite  to  austenite  and  then  to  martensite 
has  occurred.  Thus  in  these  cases,  the  hardness  is  modified  by  solid-state  phase  transforma¬ 
tion  at  locations  within  the  sample  at  which  the  temperature  reaches  the  austenization  tem¬ 
perature  (approximately  1000  K).  The  subsequent  cooling  rates  are  suflSciently  high  in  all 
cases  to  promote  the  austenite-martensite  transition.  Note  that  at  4500  mm/min,  the  depth 
at  which  hardening  occurs  is  not  measurable. 


Figure  15.  Illustration  of  microhardness  inden¬ 
tions  (100  fjm  spacing  horizontal  and  vertical 
with  50  pm  offset). 
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figure  16.  Microhardness  contour  for 
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Figure  17.  Microhardness  contour  for  4500  mm/min  beam  translation  rate. 
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Figure  18  is  a  composite  of  several  contour  plots  of  measured  microhardness  as  a  function  of 
translation  rate  for  the  beam  index  of  2.5  mm,  illustrating  that  case  depth  increases  with 
decreasing  translation  rate. 


Hardness  Sections  for  Laser  Processed  Surface 
Substrate  4340  steel,  Power=500  W,  No  Precursor 
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Figure  18.  Microhardness  map  illustrating  dependence  on  beam  translation  rate. 
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Based  on  the  cooling  rates  and  temperatures  expected  throughout  the  heat  treated  region, 
the  presence  of  a  uniform  martensitic  layer  of  some  depth  (approaching  zero  depth  as  trans¬ 
lation  velocity  increases)  was  expected.  However,  the  hardness  measurements  indicated  the 
additional  effect  of  a  soft  region  near  the  surfece  even  though  temperature  calculations  based 
on  equilibrium  heating  and  cooling  curves  for  the  ferrite-austenite-martensite  transitions 
predict  a  purely  martensitic  region.  X-ray  analysis  for  austenite  also  shows  that  there  is  no 
appreciable  austenite  in  the  layer.  The  rate  of  carbon  diflusion  in  pre-austenite  phases  has 
been  shown  to  be  of  importance  in  determining  martensite  hardness  in  cases  of  low  heating/ 
cooling  rate  laser  processing  [9],  and  this  may  be  the  case  here.  However,  there  is  also  the 
possibility  of  incomplete  austenitization  based  upon  the  high  heating  rates  for  the  laser  pro¬ 
cess  and  the  additional  effect  of  grain  growth  during  the  longer  period  of  high  temperature, 
which  is  experienced  near  the  surface.  This  is  being  investigated  further. 

Heat  Treated  4340  Processine  Trials 

In  view  of  the  above  results,  additional  studies  were  imdertaken  that  included  heat-treated 
4340  steel.  Two  sets  of 4340  steel  were  prepared.  One  set  was  heat  treated  by  bringing  to 
830°C  and  quenched  in  oil.  The  second  set  remained  as  received  in  the  non-heat  treated 
condition.  The  base  hardness  of  the  heat  treated  sample  was  680  Knoops,  and  the  base 
hardness  of  the  non-heat  treated  set  was  340  Knoops.  A  series  of  laser  tracks  were  then  run 
on  each  sample  using  the  conditions  in  Table  I.  Measurements  were  taken  of  the  melt  zone 
depth,  the  melt  zone  with  HAZ  depth,  and  the  melt  zone,  HAZ  and  tempered  zone  depths. 
These  are  shown  plotted  in  Figure  19a  and  b.  It  can  be  seen  that  the  melt  zone  depths  are 
similar  for  both  the  as-received  and  heated  treated  samples  and  they  follow  a  linear  relation¬ 
ship  with  power/velocity.  The  melt  zone  with  HAZ  depths  are  also  similar  for  both  samples 
sets,  but  they  do  not  have  a  linear  relationship  to  power/velocity.  Instead  their  linear  rela¬ 
tionship  is  to  [power/velocity]  The  tempering  zone  depth  follows  the  same  relationship 
as  the  HAZ.  This  implies  that  during  laser  surface  alloying,  a  process  which  is  generally 
confined  to  the  surface  and  near-surface  region  of  a  material,  considerations  of  heat  flow  are 
complicated  by  a  transition  from  two-dimensional  surface  effects  to  three-dimensional  ef¬ 
fects  at  depth.  Two  parallel  sets  of  microhardness  measurements  were  taken  through  the 
cross-section  of  each  track  at  50  micron  increments  with  an  automated  microhardness  tester. 
Typical  results  are  shown  in  Figure  20a  and  b,  for  a  power  of  1500  watts  at  three  speeds. 

The  heat-affected  zones  are  obvious  in  all  of  the  samples,  however,  there  are  several  differ¬ 
ences  in  the  two  sets  of  plots:  (1)  the  HAZ  in  the  as-received  samples  are  relatively  flat 
whereas  the  HAZ  in  the  heat-treated  sanples  have  distinct  slope;  (2)  The  heat-treated  samples 
display  a  tempered  region  immediately  below  the  HAZ  while  the  as-reeeived  samples  do  not. 
The  change  in  slope  of  the  HAZ  in  the  heat-treated  samples  represents  an  additional  phe¬ 
nomena  which  occurs  when  a  steel  such  as  4340  is  reprocessed  (e.g.,  essentially  laser  heat 
treated  over  an  area  that  has  been  previously  heat  treated  and  quenched).  The  region  near 
the  surface  has  ejqjerienced  two  temperature  cycles  that  reached  above  the  austenite  trans¬ 
formation  temperature  and  dropped  rapidly  to  produce  martensite.  The  microhardness  plots 
in  Fig.  20a  and  b  show  the  effects  of  a  second  cycle.  The  reduction  in  hardness  towards  the 
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Layer  +  HAZ  Depth  Squared 


surface  after  the  second  cycle  is  thought  to  be  due  to  either  grain  growth  or  a  depletion  in 
carbon  since  the  longer  times  at  temperature  encourage  difiusion. 

One  aim  of  this  segment  was  to  combine  temperature  and  time  predictions  with  the 
microhardness  data  to  estabhsh  a  rate  law  that  could  be  used  in  predicting  the  reaction  rates 
in  these  regions,  and  therefore  properties  for  a  given  set  of  conditions.  Since  transformation 
mechanisms  taking  place,  particularly  in  the  tempered  region  are  not  well  xmderstood,  it  was 
not  possible  to  develop  a  theoretical  set  of  formulations,  therefore  an  empirical  approach 
was  adopted. 
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Figure  19a.  Melt  zone  depth  as 
a  function  of power/velocity  for 
the  as-received  and  heat- 
treated  samples. 


Figure  19b.  The  square  of  melt 
zone  with  HAZ  depth  as  a 
function  of power/velocity  for 
the  as-received  and  heat- 
treated  samples. 
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Microhardness  (Knoops)  Microhardness  (Knoops) 


Figure  20a.  Microhardness 
measurements  of  as-received 
samples  processed  at  1500 
watts  and  three  speeds. 


Depth  (Meters) 


Since  tempering  has  a  time-temperature  relationship,  these  parameters 
a  rate  equation  such  as, 


Figure  20b.  Microhardness 
measurements  of  heat  treated 
samples  processed  at  1500 
watts  and  three  speeds. 


can  be  related  through 


1/t  =  A  exp  (-Q/RT) 

where  t  is  time,  A  is  a  dimensionless  constant,  R  is  the  gas  constant,  Q  is  an  activation  energy 
representing  the  process,  and  T  is  temperature.  The  effectiveness  of  tempering  is  often 
expressed  in  terms  of  hardness  and  so  this  equation  can  be  rewritten  as  [8] 

Hardness  =  f[tA  exp  (-Q/RT)] 

The  simplicity  of  this  equation  is  dependent  upon  certain  assumptions  such  as:  A,  which 
reflects  the  microstructure  is  independent  of  tenq)erature  and  Q  is  independent  of  hardness. 
Hence,  the  activation  energy  Q  for  a  reaction  can  be  obtained  from  data  that  is  gathered  as  a 
fimetion  of  temperature.  Some  typical  activation  energies  of  interest  are  given  in  Table  V. 
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TABLE  V 


REACTION 

ACTIVATION  ENERGY  (CAL/MOLE) 

Carbon  diffusion  in  fee  iron 

34,000 

Carbon  diffusion  in  bcc  iron 

29,300 

Iron  diffusion  in  fee  iron 

64,000 

Iron  diffusion  in  bcc  iron 

57,500 

Cementite  growth  at  300“C 

48,500-j/moIe  jl-6096 

The  larger  the  activation  energy,  the  slower  the  reaction. 

Measurement  of  the  activation  energy  Q  in  the  tempered  region  was  imdertaken  using  the 
sets  of  microhardness  data  for  the  two  base  materials.  In  order  to  represent  the  time-tem¬ 
perature  relationship,  the  depths  of  the  molten  zone  and  HAZ  were  fit  to  an  existing  empiri¬ 
cal  equation  [9].  Using  these  corrections  on  the  model  material  parameters,  the  time  period 
over  which  each  microhardness  mark  location  was  at  an  elevated  temperature  was  calcu¬ 
lated.  The  inverse  of  this  time  was  then  plotted  logarithmically  versus  the  change  in 
microhardness  with  time  at  selected  peak  temperatures.  Figure  21  shows  this  plot.  An 
activation  energy  of  3  5  kcal/mol  can  then  be  obtained  from  Figure  2 1 .  This  is  in  line  with  the 
diffiision  of  carbon  in  iron,  and  although  it  is  based  on  an  analytical  calculations  gives  an 
indication  of  the  controlling  phenomena.  It  also  indicates  the  value  and  the  need  for  an 
accurate  computational  model  that  is  validated  by  experimental  data.  Such  an  effort  was 
undertaken  in  the  third  phase  of  this  project. 


Figure  21.  The  slope  of  hardness/time  as  a  function  of  the  ln(l/T)  where  T  is  temperature. 
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ni.  F.xnerimental  and  Computational  Thermal  Studies 


Thermocouple  measurements  were  used  to  determine  the  heating  and  cooling  rates  during 
laser  surface  processing.  Because  mechanical  measurements  were  being  made  on  a  micro¬ 
scale,  and  regions  of  interest  may  span  only  150  microns,  temperature  measurements  on  the 
same  scale  were  desired.  One  basic  objective  of  the  measurements  was  to  determine  directly 
the  thermal  cycle  in  the  near  melt  region.  In  order  to  accomplish  this,  thermocouples  must 
be  placed  in  this  region  in  such  a  way  that  they  represent  the  actual  temperature  in  this 
region.  This  was  accomplished  for  the  short  span  of  time  before  the  thermocouple  reaches 
its  melting  temperature.  Figure  22  is  a  comparison  of  temperature  measurements  made  at 
the  same  sample  location  for  trials  of  differing  index.  The  temperatures  for  the  small  index 
case  are  noticeably  higher  than  the  large  index  case.  In  some  cases  the  temperature  may  be 
sufficient  to  provide  an  annealing  process  for  previously  processed  material.  This  has  been 
observed  in  some  samples  with  very  small  beam  index  (.25  mm). 


Figure  22.  Comparison  of  thermocouple  measurements  for  3  mm  and  2  mm  beam  index  values, 
(a)  small  index  case,  and  (b)  large  index  case. 
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Several  computational  models  were  developed  for  predicting  the  three-dimensional,  time- 
dependent  temperature  in  the  processed  samples  and  were  used  to  investigate  the  influence 
of  various  operating  conditions  on  the  process.  One  was  a  fimte  volume  model  that  can 
accurately  predict  multi-species  convective  transport,  and  the  other  was  a  fimte  element 
model  that  can  be  used  to  predict  thermally  induced  stresses.  Although  there  was  some 
overlap  between  the  models,  they  provided  different  physical  information,  and  the  objective 
was  to  calculate  the  temperature  and  species  transport  with  the  finite  volume  model,  and  use 
this  as  input  to  the  finite  element  stress  model. 

The  ultimate  goal  of  the  model  development  was  to  determine  the  most  important  character¬ 
istics  of  the  laser  surface  processing;  melt  pool  dimensions  and  chemical  composition  distri¬ 
bution  being  two  primary  examples.  However,  the  importance  of  the  time  dependence  of 
temperature  in  determining  sohd  phase  formation  drove  the  initial  efforts  at  modeling  the 
process  towards  solving  the  simple  three-dimensional  heat  equation,  i.e.,  the  time  dependent 
difiusion  equation.  While  the  actual  temperature  field  can  be  determined  only  by  solving  the 
full  set  of  coupled  conservation  laws  and  accurately  assessing  the  physics  of  the  phase  trans¬ 
formation  at  any  phase  transition  boundaries,  the  heat  conduction  equation  is  useful  as  a  first 
approximation  to  the  temperature  field  in  the  sample. 

The  thermal  model  is  an  aid  in  predicting  cooling  rates  for  a  material  that  has  been  laser 
processed.  Since  there  are  several  process  variables  that  must  be  considered  (beam  power, 
beam  speed,  overlap  index),  a  model  can  be  used  to  determine  cooling  rates  for  a  given  set  of 
conditions.  In  any  given  track,  in  addition  to  the  phase  transformations  within  the  melt  pool, 
there  can  be  significant  heating  in  the  heat  affected  zone  (HAZ),  i.e.,  the  regions  adjacent  to 
the  melt  pool  that  are  heated  to  temperatures  close  to,  but  below  the  melting  point.  Since 
there  is  a  solid  state  phase  transition  at  1000  K  in  steel  and  melting  typically  occurs  aroimd 
1700  K,  it  is  of  interest  to  know  which  regions  in  the  vicinity  of  the  melt  pool  are  heated 
above  the  transition  temperature.  The  model  can  determine  this  region.  The  model  can  also 
be  used  to  determine  the  cooling  rates  of  the  HAZ,  in  order  to  estimate  the  phase  distribution 
around  the  melt  pool. 

Since  the  surfece  processing  utilizes  a  cyclic  process,  there  are  multiple  beam  tracks  offset 
by  a  given  index,  and  it  is  possible  that  the  HAZ  of  one  track  may  overlap  the  melt  pool  of  a 
prior  track  and  yield  an  additional  soM  state  phase  formation.  It  is  important  to  determine  if, 
for  a  given  set  of  conditions,  there  is  enough  overlap  to  heat  processed  regions  to  tempera¬ 
tures  above  the  lower  transition  temperature,  and  hence  possibly  transforming  it  into  a  dif¬ 
ferent  phase  than  that  immediately  following  the  original  transition  fi-om  liquid  to  solid. 

For  the  finite  volume  model,  the  commercial  computational  fluid  dynamics  code,  CFDACE+ 
(CFD  Research  Corporation,  Huntsville,  AL),  was  used  to  solve  the  energy  equation.  The 
model  is  run  on  a  PC  using  the  W^idows  NT  (v4.0)  operating  system.  CFD  ACE  has  the 
advantage  that,  in  the  fixture  when  momentum  and  mass  transport  will  be  included,  it  can 
solve  the  coupled  equations  associated  with  fluid  flow,  and  hence  can  be  used  to  model  the 
miying  and  distribution  of  species  within  the  melt  pool.  In  the  current  model  there  is  no 
consideration  of  mass,  momentum,  or  species  transport. 
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Figure  23  depicts  the  domain  of  the  problem.  The  size  of  the  3D  computational  domain  is 
limited  by  memory.  The  scale  of  the  laser  requires  a  fine  grid  in  the  vicinity  of  the  beam  in 
order  to  model  the  beam  profile  and  steep  thermal  gradients  around  the  beam.  Because  the 
model  is  memory  limited,  each  computational  grid  is  tailored  to  a  specific  set  of  processing 
conditions.  The  grid  in  Figure  23,  for  example,  is  designed  to  handle  three  tracks  with  an 
overlap  index  that  is  1/2  the  width  of  the  beam.  If  the  five  densely  packed  regions  in  the 
center  of  the  grid  are  hypothetically  numbered  1  through  5  in  order,  then  the  first  track 
would  be  centered  on  2  with  the  edges  of  the  beam  located  in  1  and  3.  The  second  track 
would  similarly  be  centered  on  3,  with  the  beam  edges  located  in  2  and  4.  Other  processing 
parameters  are  used  to  define  other  unique  grids. 


Fi^re  23.  Computational  Grid  (251  points  in  the  direction  of  the  beam  movement,  16 
points  vertical,  and  62  points  perpendicular  to  beam  movement). 

The  grid  spacing  is  varied  in  both  the  y  direction  and  z  direction  in  order  to  pack  points  near 
the  beam.  The  x  direction  cell  spacing  is  100  pm.  the>’  direction  cell  spacing  varies  from 
2000  pm  at  the  edge  of  the  model  to  a  minimum  of  150  pm  in  the  densely  packed  "track" 
regions.  The  z  direction  cell  spacing  varies  from  2000  pm  on  the  side  away  from  the  beam 
(adjacent  to  the  heat  sink)  to  30  pm  at  the  processing  surface.  Memory  constraints  limit  the 
model  to  a  volume  of  2.5  cm  (x,  I  =  251)  x  2.5  cm  (y,  J  =  62)  x  6  mm  {z,  K  =  16). 

The  boundary  conditions  are  selected  to  closely  approximate  actual  processing  conditions, 
i.e.,  convective/radiative  and  or  fixed  temperature.  Since  a  shield  gas  is  typically  used,  the 
convective  heat  transfer  coefiicient  is  somewhat  higher  than  “natural,”  30  W/m^  K,  to  ac¬ 
count  for  the  forced  flow  of  gas  around  the  sample.  All  five  exposed  sides  use  this  value.  The 
sample  is  mounted  to  a  water-cooled  heat  sink,  so  the  bottom  side  could  be  set  to  a  fixed 
temperature.  However,  experience  shows  that  the  water  flow  can  only  remove  some  of  the 
energy  building  into  the  heat  sink,  and  thus  an  alternative  boundary  condition,  is  to  use 
convection  (no  radiation)  with  a  high  (500  W/m^)  heat  transfer  coefiicient. 
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Two  beam  profiles  that  were  modeled;  an  elliptic  Guassian  and  rectangular  tophat,  both  with 
a  possible  temperature  dependent  absorption  coefficient. 

The  Guassian  beam  is  of  the  form: 


while  the  rectangular  tophat  is  of  the  form: 


y-ydf) 

.  yr  ) 


qb(x,y,t)  =  alg  when  \x  -  x^it)  <x^ax\A\y -yg{t)\<y^ 

=  0  otherwise 


where  xjii)  and>’g(0  are  the  coordinates  of  the  center  of  the  moving  beam.  For  the  Guassian 
beam  profile,  x^  and>/^  are  the  beam  radii  to  the  1/e  point  in  their  respective  directions  for  the 
rectangular  tophat  profile,  x^  is  1/2  the  beam  width  in  the  x  direction  and  y^  is  1/2  the  beam 
width  in  the  y  direction. 


is  the  beam  intensity  at  the  center  of  the  beam  and  can  be  found  fi-om  the  beam  power,  P, 
using  p 

^0  ~cx  y 

r-y  r 

where  c  is  7t  (3.14159)  for  the  Guassian  profile  and  4  for  the  rectangular  tophat  profile. 

Solving  the  time  dependent  heat  equation  requires  only  the  density,  thermal  conductivity, 
and  specific  heat.  The  density  is  assumed  constant,  while  the  conductivity  and  specific  heat 
are  fimctions  of  temperature.  Additionally,  since  the  model  does  not  incorporate  phase  changes 
ejqjlicitly,  the  latent  heats  of  formation  are  included  in  the  specific  heat.  Figure  24  illustrates 
the  entMpy  distribution  for  the  shown  temperature  profile  using  the  specific  heat  of  Table 
VI. 


Figure  25  shows  the  isotherms  of  the  temperature  field  at  various  times  during  processing 
with  a  low  power  (200  W)  beam.  Note  that  even  in  this  case  of  minimal  energy  input,  there 
is  evident  preheating  of  the  sample  for  subsequent  laser  tracks.  This  yields  an  increasing 
maYimiim  temperature,  and,  together  with  the  thermal  boundary  conditions,  an  asymmetry 
in  the  temperature  distribution. 

Figure  26  compares  the  measured  temperature  profile  during  processing  using  a  moderate 
power  (500  W)  beam  with  modeled  values  using  a  rectangular  tophat  beam  profile  with  the 
same  power.  Calculated  heating  and  cooling  rates  match  measured  values  very  well.  Slight 
differences  can  be  attributed  to  measurement  technique. 

In  addition  to  the  volume  control  model,  a  finite  element  model  was  constructed  using  the 
software  package  COSMOS/M.  The  objective  for  this  model  is  to  predict  (or  input  from  the 
CV  model)  the  temperature  field  and  resulting  stresses  during  and  after  LISI  processing. 
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Typical  Enthalpy-Temperature  Curve 


Figure  24.  Enthalpy  as  a  function  of  temperature. 


TABLE  VI 

Typical  Specific  Profile  for  a  "Steel" 


Tenperature  Range  (K) 

Specific  Heat  (J/kg  K) 

1-1700 

460 

1700-1900 

1100 

1900-3000 

490 

3000-3200 

24000 

3200- 

IL-6096 

30 


Figure  25.  Coni. 

Calculated  temperature  contours. 

(d)  end  of  2nd  track, 

(e)  beginning  of  3rd  track, 

(f)  end  of  3rd  track. 


Figure  26.  Comparison  between  modeled  and  measured  time  dependent  surface  tempera¬ 
tures  at  two  different  locations.  The  solid  lines  represent  thermocouple  measurements. 

The  dashed  lines  represent  modeled  temperatures  at  the  same  location.  The  dotted  lines 
represent  modeled  values  300  microns  (the  diameter  of  the  thermocouple  bead)  closer  to 
and further  from  the  beam. 
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The  problem  to  be  modeled  was  similar  to  the  sample  coupons  used  in  the  processing  trials, 
with  dimensions  3  cm  x  3  cm  x  0.6  cm.  The  grid  in  the  vertical  dimension  is  comprised  of  10 
elements  that  have  nonuniform  spacing  varying  from  .29  microns  at  the  top  surface  to  1764 
microns  at  the  lower  surface.  There  are  50  elements  in  each  of  the  other  directions.  Bound¬ 
ary  conditions  at  the  side  and  upper  surfaces  are  a  combination  of  convection,  h  =  30  W/m^ 
K  (ambient  temperature  =  20'’C)  and  radiation.  The  bottom  surface  is  adjacent  to  a  water 
cooled  heat  sink  and  is  modeled  by  h  =  500  W/m^  K  and  temperature  =  20‘’C.  The  distribu¬ 
tion  of  energy  in  the  laser  beam  is  modeled  as  an  elliptic  Gaussian  and  the  surface  has  an 
absorption  coeflScient  of  0.54.  Figure  27  shows  a  comparison  of  temperature  responses  at 
various  locations  at  the  top  of  the  surface.  Note  that  although  the  track-to-track  tempera¬ 
ture  rise  is  evident,  this  model  does  not  capture  the  high  heating  and  cooling  rates  correctly. 
There  is  an  oscillatory  behavior  that  is  especially  noticeable  during  the  initial  heating  of  the 
point  nearest  the  beam.  This  is  a  common  problem  with  explicit  time  solutions,  and  is  one 
reason  the  more  robust  CFD-ACE  model  with  a  fiilly  implicit  solver  was  selected  for  primary 
use  in  this  project. 


Figure  27.  Calculated  Temperature  Response. 
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Figure  28  shows  the  calculated  temperature  response  at  the  top  three  grid  points  during 
multi-track  processing.  Note  that  at  this  translation  rate,  all  of  the  points  shown  are  heated 
above  the  austenitization  temperature.  Thus  the  material  at  these  locations  will  be  trans¬ 
formed  to  martensite  and  will  exhibit  higher  hardness  than  the  substrate.  This  data  can  be 
examined  to  identify  a  hardening  depth  that  can  be  compared  to  measurements. 


Computed  Temperature-Time  Curves  for  Near-Surface  Region 


Figure  28.  Temperature  response  at  the  top  of  the  processed  surface  during  multi-track  processing. 
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SUMMARY 


In  studying  laser  surface  alloying,  there  are  two  significant  areas  of  concern:  the  alloyed 
surfece,  and  the  base  material.  The  alloying  process  itself  presents  myriad  potential  phenom¬ 
ena  depending  upon  the  elements  and  compoimds  used  and  their  afiSnity  for  chemically  re¬ 
acting  with  the  base  material  or  each  other.  The  base  material  or  region  below  the  surfece 
alloy,  presents  another  set  of  interesting  phenomena  in  that  it  can  experience  numerous  met¬ 
allurgical  events  such  as  martensitic  transformation,  tempering,  grain  growth,  carbide  growth 
etc.  The  most  desirable  approach  to  identifying  the  total  resultant  properties  of  such  a  laser 
surface  alloyed  material  would  be  to  exercise  a  computational  model.  This  study  attempted 
to  begin  the  development  of  such  a  model  while  performing  concurrent  experiments  to  iden¬ 
tify  the  significant  chemical  reactions  and  metallurgical  phenomena  that  must  be  considered. 

A  systematic  approach  was  adopted  to  investigating  the  addition  of  materials  into  the  mate¬ 
rial  surfece.  When  considering  steel,  it  was  determined  that  the  effects  of  elemental  addi¬ 
tions  such  as  nickel  and  chromium  can  be  predicted  using  existing  phase  information  for  the 
stainless  steels.  However,  when  compoimds  are  included  in  the  additive,  chemical  reactions 
occur  that  create  important  deviations  from  the  alloy  composition.  These  must  be  consid¬ 
ered  individually  since  the  reactions  are  dependant  both  upon  the  relative  amount  of  the 
additive  and  the  particular  additive  itself  A  good  example  from  this  study  is  the  result  of 
adding  WC  to  steel.  Decomposition  of  the  WC  into  other  compounds  of  tungsten  and  the 
addition  of  carbon  to  the  matrix,  produces  a  unique  set  of  properties  within  the  surface 
material.  Although  the  properties  such  as  hardness  and  wear  resistance  are  improved  by  the 
WC  addition,  the  microstructural  path  to  obtaining  those  properties  must  still  be  understood 
and  predictable. 

In  steels  in  particular,  the  region  below  the  alloyed  surface  displays  a  variety  of  mechanical 
properties  which  depend  upon  the  processing  parameters.  The  energy  input  from  the  laser 
determines  the  time  and  temperatures  which  are  experienced  by  the  steel  and  therefore  influ¬ 
ences  the  solid  state  transformations.  Although  it  has  been  thought  that  the  speed  of  the 
laser  mitigates  many  of  these  transformations,  this  is  not  the  case.  The  transformation  to 
austenite  and  then  to  martensite  occxirs  within  the  subsecond  time  frame  of  the  laser  pass, 
and  more  interestingly,  the  region  beneath  the  heat  affected  zone  also  has  sufficient  time  to 
temper.  These  are  effects  which  have  been  found  as  a  result  of  this  investigation  and  bear 
fiirther  study  to  distinguish  between  the  influence  of  the  temporal  cycle  on  grain  growth, 
carbon  difiusion,  carbide  formation,  and  phase  transformation.  Processing  parameters  that 
can  both  produce  a  layer  with  selected  properties  while  at  the  same  time  encouraging  benefi¬ 
cial  phase  transformations  below  the  layer  exist,  but  within  a  narrow  window  and  one  that 
needs  to  be  clearly  examined  and  identified. 

Utilization  of  a  computational  model  for  the  prediction  of  layer  and  sub-layer  phases  is  a  far 
more  desirable  approach  than  running  a  multitude  of  laboratory  experiments.  However,  as 
the  experimental  studies  elucidated,  the  reactions  and  transformations  which  occur  during 
and  subsequent  to  processing  are  con^licated  and  varied.  This  study  began  attacking  that 
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large  problem  with  a  3D  model  that  eould  produce  the  melt  pool  and  the  cooling  rates.  Since 
the  surface  processing  itself  is  cyclic,  the  indexing  of  the  tracks  was  also  considered  since 
one  track  overlaps  the  solidification  region  and  HAZ  of  the  adjacent  track.  The  grid  spacing 
was  varied  in  two  directions  in  order  to  correctly  model  the  intensity  of  the  small  beam. 
Calculated  temperature  contours  for  repetitive  tracks  were  ft)und  to  closely  match  data  taken 
fi'om  high  speed  temperature  measurements.  Therefore  the  regions  which  exceeded  both  the 
liquid  and  solid  transformation  temperatures  was  identified  for  an  actual  experimental  sample 
configuration.  The  influence  of  a  previous  laser  track  upon  the  subsequent  track,  was  also 
distinguished.  This  is  a  significant  first  step  in  understanding  and  predicting  the  final  proper¬ 
ties  of  a  laser  surface  alloyed  material. 
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